function img=img_polar(I,radius,angle)


   if nargin < 1
      error('Porfavor especifica una imagen');
   end
   
   
   [rows,cols] = size(I);
   cy          = round(rows/2);
   cx          = round(cols/2);
   
   if exist('radius','var') == 0
      radius = min(round(rows/2),round(cols/2))-1;
   end
   
   if exist('angle','var') == 0
      angle = 360;
   end
  
   img1 =zeros(rows,cols);
   i     = 1;
   
   for r=0:radius
      j = 1;
      for a=0:2*pi/angle: (2*pi) - (2*pi/angle)
         %pcimg(i,j) = img(cy+round(r*sin(a)),cx+round(r*cos(a)));
         if(r>50 && r <100)
           %img1(cy+ceil(r*sin(a)),cx+ceil(r*cos(a)))=1;
           img1(cy+round(r*sin(a)),cx+round(r*cos(a)))= 1;
           
        
           img1(cy+round(r*sin(a))-1,cx+round(r*cos(a))-1)= 1;
           img1(cy+round(r*sin(a))-1,cx+round(r*cos(a)))= 1;
           img1(cy+round(r*sin(a))-1,cx+round(r*cos(a))+1)= 1;
           img1(cy+round(r*sin(a)),cx+round(r*cos(a))-1)= 1;
           img1(cy+round(r*sin(a)),cx+round(r*cos(a))+1)= 1;
           img1(cy+round(r*sin(a)+1),cx+round(r*cos(a))-1)= 1;
           img1(cy+round(r*sin(a)+1),cx+round(r*cos(a)))= 1;
           img1(cy+round(r*sin(a)+1),cx+round(r*cos(a))+1)= 1;
           
         end
         
         j = j + 1;
      end
      i = i + 1;
   end
   img=img1;
   
end